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Abstract 

A discrete rotation algorithm can be apprehended as a parametric map f a from Z[i] to 
Z[i], whose resulting permutation “looks like” the map induced by an Euclidean rotation. For 
this kind of algorithm, to be incremental means to compute successively all the intermediate 
rotated copies of an image for angles in-between 0 and a destination angle. The discretized 
rotation consists in the composition of an Euclidean rotation with a discretization; the aim of 
this article is to describe an algorithm which computes incrementally a discretized rotation. The 
suggested method uses only integer arithmetic and does not compute any sine nor any cosine. 
More precisely, its design relies on the analysis of the discretized rotation as a step function: 
the precise description of the discontinuities turns to be the key ingredient that will make 
the resulting procedure optimally fast and exact. A complete description of the incremental 
rotation process is provided, also this result may be useful in the specification of a consistent 
set of definitions for discrete geometry. 

keywords : Discrete Rotations, Discrete Geometry , Computational Geometry. 


1 Introduction 

The translation of the fundamental concepts of the Euclidean geometry into 7L n comprises the field 
of discrete geometry. As this theory of geometry is particularly suitable for combinatorial images 
and other data manipulated by computers IKR04I , it would be interesting to provide a set of efficient 
algorithm for this theory that uses onlv integer-arithmetic: as this was suggested in (Rev nj. 

Several attempts have been realized by various authors that wished to deliver back some prop¬ 
erties of the Euclidean rotation to the discretized rotations widely used in computer graphics. A 
review of various attempts may be found in IAnd92j . 

In this paper, we present an algorithm which is incremental: it successively computes all the 
rotated images according to the an increasing sequence of angles (starting from 0 to 27 t). Notice 
that the set of rotated images is finite on a finite picture, this allows practically to compute all the 
intermediate rotated images. Moreover, the suggested procedure is sound and accurate: it returns 
exactly the same results as the result provided by the discretized for the same angle. The procedure 
does not use any sine nor any cosine, thus there is no precision problem due to the floating point 
arithmetic. Also, the algorithm is fast: to compute incremental rotations the algorithm computes 
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only 0(m 3 — log(m)) operations, instead of 0(m 5 ) as need the naive algorithm. For incremental 
rotations, the complexity of this algorithm if it uses pre-calculated tables, 0(m 3 ) is optimal: The 
algorithm updates only the necessary pixels and only consider the necessary angles. Finally, due 
to the fact the algorithm uses configurations that can be stored with very few states on the plane, 
we believe it is a good candidate for parallelization. 

After a brief review of the motivations, and after the essential preliminary definitions, we proceed 
to a characterization of the discontinuities of the rotation process. Indeed, we will explain how to 
code the angles where the discontinuities happen. Also, with integer arithmetic only, we will learn to 
perform essential operations on the encoded angles. Naturally, a few technical lemmas are required 
to set up all this framework properly. Once this has been set, we will analyze the alterations that 
occur in the configuration at the discontinuities. Strengthened by previous results, we will then be 
ready to build the incremental discretized rotation procedure. The last section will be dedicated 
to various extensions and miscellaneous details related to the theory that may lead to a better 
understanding of the discretized rotation process. 

2 Groundwork 

The first sections present the fundamental ideas, definitions and lemmas that matter to fully un¬ 
derstand the algorithm. In this section, we review the motivations and precise some vocabulary. 

2.1 Motivations 

The history of discrete geometry begins with the common will to give birth to an algorithmic the¬ 
ory of the geometry in the discrete spaces that would be consistent with the Euclidean geometry. 
We believe that a unified theory would provide a better understanding of both universe. Discrete 
rotations comprise the famous examples that have strengthened the idea that discrete and contin¬ 
uous spaces may be radically different. A review of the differences can be found in the prologue of 

EE053. 

More recently, discretized rotations have been an important issue in water-marking community. 
Water-marking algorithms that were robust under rotation were sought for by various teams the 
community; and it supports many discussions. More generally, the problem of finding algorithms for 
classification and recognition of patterns that are robust under discrete rotations is still a not-trivial 
issue in the conception of pattern related algorithms. The algorithm presented use similar principle 
has the one suggested by fABC + 04| ; the basic idea is to compute a table of the discontinuities of the 
discrete rotation process, and to use this table for the rotation. Our paper provides a more accurate 
description of the discontinuities, as-well as a way to compute the rotation using integer-arithmetics 
only. 

The algorithm was designed to comprise the field of discrete geometry: our main problem is to 
provide a comprehensive and simple definition of algorithm for the discrete rotation, independantly 
of any specific usage, and that can be computed efficiently. 

The usual discretized rotation algorithm is a bit chaotic when iterated: This is well illustrated 
by the following story: It has been implemented a physical simulation of the solar system which 
for specific reasons was coded with fixed precision. When the system has been tried, its inventors 
were surprised that it first conclusion was to predict a fatal collision in-between the earth and the 
sun in less than ten years; naturally there was a “bug”, and it was imputable to composition of 
rounding errors during the rotation process. More generally, decades of computer programmers 
have known by experience that they should avoid to compose rotations: The accumulation of the 
resulting errors may produce an unwelcome result. Some aspects of the dynamical system that is 






formed by the iterated action of a discretized rotation have already been studied; see |Pog03| for 
iteration of 7r/4-rotations. 

2.2 Conventions 

We work in the complex plane C, where Z[i\ denotes the set of Gaussian integers. Let m be a 
positive integer. We denote by Z[z] \ m the set of Gaussian integers whose module is at most m; 
Z[i]\ m = {z E Z[z], \z\ < m}. Real and imaginary parts are denoted 3i(z) and 3f(z). Let x be a real 
number. We recall that the floor function x 1 —> |_®J is defined as the greatest integer less or equal 
to x. The rounding function is defined as : [x] = [x + 0.5J ; we also define the map x 1 —> {x} by 
{x} = x — [x]. These maps are extended to complex numbers, by applying them independently on 
the real part and on the imaginary part. Let TL be set of complex number that have a semi-integer 
coordinate, in the general case Ti denotes the set of discontinuity points of the operator x 1 —> [x] ( 
in other terms Ti = (M x (Z + {^})) U ((Z + {^}) x M) ). 

Let a denote a angle, i.e. an element of A = M/(27tZ). The Euclidean rotation r a is the 
bijective isometry of C: z 1 —> ze ia . The discretized rotation [r a ] is precisely defined as the successive 
computation of the Euclidean rotation of angle a and of the discretization operator z >—> [z\. Thus, 
for each z of C, [r a (z)] = [ze ia ]. Remark that |[r a (^)]| < \z\ + and that | \r a (z)\ — r a {z)\ < 

A configuration is a mapping from Z[i] to Z[i\. Let m be a positive integer. A configuration C 

such that for all z E C, lim^^oo —-( e lOL z ) —> 0 is called an a-rotation map. Thus, given a 

real a, the discrete rotation [r Q ] is a configuration. A partial configuration of radius m is a mapping 
from Z[i]| m to Z[i]. Each configuration induces a partial configuration. In this paper, we work on 
on partial configurations (denoted by [r a ]| m ) induced by discretized rotations. Precisely, we study 
the mapping p m : a *—> [r Q ]| m . Since Z[i\\ m is finite and |[r a (z)]| < \z\ + 1, the set {p m (a),a E R} 
is hnite. 

In this paper, we produce an exact incremental rotation algorithm which, given an integer m 
successively produces all the values p m ,(a) for a E [0, 27t], in the order where they are reached (mov¬ 
ing from 0 to 2tt). Informally, it is a ” video” algorithm which exhibits the successive configurations 
obtained along the rotation. 

3 Hinge Angles 

3.1 Definitions 

The principal element that have influenced the design of the algorithm is a precise study the function 
p m . We shall prove that p rn is a stair function(piecewise constant function). Thus, a precise study 
of its discontinuity steps gives the ability recover the whole function. The discontinuities also 
correspond to the only updates of the configuration stored in the algorithm we will describe. 

Definition 1 An angle a is a hinge angle if there exists a source point z s in Z[i] such that the 
destination point z r i = z s e ia has a (proper) semi-integer component (i.e. z&^TL) . For each hinge 
angle, the source point and the destination point form a generating pair. 

On a topological point of view, the hinge angles are the discontinuity points of the map a e-> [r a ]. 
It is immediate from the definition that an angle is a hinge angle if and only if there exists 
integers p, q, k such that 2 q cos (a) + 2psin(a) = 2 k + 1. 

The value |z s (a)| 2 is called the order of a. We can note that if (z s , Zd) generates a then (z s , Zd) 
generates —a.Note, also that if ( z s ,Zd ) generates a then (iQz s ,i®Zd) generates the same angle, for 
any Q E {0,1, 2,3}. We also define S a the set of^source points of a: S a = {z E Z[i]\e ia z E Ti} 






Figure 1: An hinge angle and its generating pair for a = a(9, —4, 7). 


3.2 Fundamental Lemmas 


The Pythagorean angles , as seen in XHO-ll or |Vos93j . are such that a = arctan(a/6) where a and 
b are issued from a Pythagorean triple (a,b,c) £ N 3 (such that a 2 + b 2 = c 2 ). c will be called the 
radius of the Pythagorean angle. An angle is Pythagorean if and only if its cosine and sine are 
rational. 


Lemma 1 Hinge angles and Pythagorean angles form disjoint sets. 

Proof: Assume that there exists an angle which is both Pythagorean and hinge. By definition, 
as it is a hinge angle, there exists a Gaussian integer z s of Z[i] which is transformed by rotation 
in a point of hi. We can easily check that: z s )} = More exactly, iR(zd) = 

5P(z s ) cos(cc) — sin(o:). If a is Pythagorean, cos(a ) = | and sin(a) = -, where ( a,b,c ) is a 

primary Pythagorean triple: i.e. belongs to Z and a 2 + b 2 = c 2 , gcd(a,b,c) = 1. {K(zrf)} can be 
written as ^ where n as an integer, while c is odd. But this is contradicting: = | □ 

A corollary of this lemma is that for a hinge angle, at least one of cos(a) and sin(a) is an 
irrational quadratic number. 

Let p,q,k be a triple of integers, such that p 2 + q 2 > \k + ^|, we state r 2 = p 2 + q 2 , and 
A = \Jr 2 — (k + ?j) 2 we define the angle a(p, q, k) by the equality e ia ( p,q,k \p + qi) = k + | + Xi. 
Thus for this angle ((p + qi), (k + ^ + Xi)) forms a generating pair. 

Lemma 2 (Coding of Hinge Angles) Let ( p , q, k) and (p', q ', k’) be two distinct generating triples 
such that a(p, q, k) = a(p', q', k’). Then det ^, C \ = 0. 


Proof: The proof is straightforward: Let a = a(p,q,k), thus 2psin(a) + 2gcos(a) = 2k + 1, 
and similarly; 2p' sin(a) + 2q' cos(a) = 2k! + 1. Consider these two equations as a linear system 


of cos(a) and sin(a). If det( 


p q 
p' q' 


) 7 ^ 0 then cos(cc) and sin(a) have to be rational, and the 


angle have to be Pythagorean. But this would contradict lemma ^1 Therefore (p, q) and ( p',q') are 
colinear. □ 

















The lemma implies that arg(z s (a)) mod (vr/4) and arg(zd(a)) mod (n/4) are uniquely defined 
for an angle a. Thus, they are important characteristics of the hinge angle. An immediate corollary 
of Lemma [3 is that all generating pairs of a(z s ,Zd) are necessarily of the form ( kz s ,kzd ) or of the 
form (ik z s ,ik Zd) with k E R. 

For a triple (p, q, k ) that generates an angle a,(p, q. k) is called primary if it is minimal p 2 + q 2 
among the generating triples of a. Obviously , there exists for each hinge angle a unique primary 
generating triple. 

Lemma 3 (Primary Generating Triple) Let (p, q, k) be the primary generating triple of an 
angle a; the set of generating triples of a is {((2 n + 1 )p, (2 n + 1 )q, (2 n + 1 )k + n), for n € Zj. 

Proof: It is evident that if z s = (p + qi) E S a is the source point issued from the primary 
generating triple (p, q, k ) then (2 n + 1 )z s E S a , since {3R(2?z + 1 )e lol z s } = \ for any n in Z. Also we 
can notice that for any n E Z, this is another generating triple ((2n + 1) p, (2 n + 1) q, (2 n + 1) k + n) 
of the same angle. It is also evident that for any n', n" E Z, (2 pn' 1 2 qn', n") cannot be a generating 
triple of a since e ia (p + qi) E Ti implies that e ia (2pn' + 2 qn'i) ^ hi. Now, we can notice (3p + j3qi 
with (3 E R \ Z cannot be a source for arithmetical reasons: /3(p + qi) ^ Z [i\. Now, assume there 
exists | E Q \ Z and k! E Z such that (|p, | qi , k!) is a generating pair of a, this would implies that 
|(p + qi) is a Gaussian integer and thus gcd(p,q) should be divisible by s (as r is not). Moreover 
$t(e ta (?-p + ^qi)) = 2fc 2 +1 . This can only happen if 2k' + 1 is also divisible by s. Now p, q, 2k + 1 
are all divisible by n, thus the generating pair is not primary, and this contradicts our hypothesis. 
Thus, f3(p + qi) ^ S a for /3 E M \ Z. Thus, consequently to previous lemma and to these points, 
there is no other generating triples for the angle a than the listed ones. □ 


3.3 Main Properties 

From the previous lemmas, we can state the following properties that are useful for the rotation 
algorithm. 

Proposition 1 Any hinge angle can be uniquely described by its primary generating triple. 

This proposition is actually a corollary of Lemma |3 and |3 

Proposition 2 The number of hinge angles of order at most m is lower than 8m 3 . 

This is easily proven by the fact that \p\ < m, \q\ < m, \ k + < m 

Upper bounds on the number of possible hinge angles can be found in lABC + fl4j . This formula 
can be slightly refined by using r 2 {k) that represents the number of decomposition of an integer 
as the sum of two squares. The upper bound on the number of possible hinge angles is then: 
#(*4w|m) < YZi r 2 {i)[Vi— < (E£i r 2 (i)) Y\/rn — < (2m) 3 . Finally note that there is twice 

more rotation maps in p m (a) than the number of hinge angles in Anlm,- 

Proposition 3 The elements of S a forms 4 rays: S a = UQ e { 0 ,i i 2,3}{(2^ + 1 )i® (p + qi),l £ N} 


This proposition is another consequence of EL 
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4 The Algorithm 


Schematically, the algorithm is structured in two main parts: During this first period of the al¬ 
gorithm, it starts with the enumeration of all the hinge angles whose order is smaller than the 
maximum order of the points of the picture. The angles are not represented by floating point val¬ 
ues but only via their associated triples (three small integer numbers). This encoding of the angles 
provides a way to recover all the necessary information required to transform one image into the 
next image during the incremental rotation process. This first part of the algorithm, can be seen 
as a process similar to the generation of a sine table and, ideally, it could be done once for all. 

The second part of the process consists in a loop through the cycle of hinge angles. This loop 
applies successively the small transformations that are required to pass from p m (a) to Pm(a') ( 
where a! is the successive representative angle ). Practically, the algorithm stores the map p m (a), 
this allows to know directly the position of the image through the transformation [r Q ]| m . However, 
for various usages, such rotations in cellular automata, the user should may prefer variant encoding 
of the transformation such as rotation configurations (maps from Z[i] to a finite set) jNRQal . The 
principle of this algorithm is translatable on these kind of configurations. 

4.1 Enumerate and Sort Hinge Angles 

To enumerate all hinge angles in An\mi it sufficient to enumerate all the Gaussian integers whose 
module smaller than m (as source point), and to consider all possible semi-integers whose absolute 
value is smaller than m. This requires a time of 0(m 3 ). Some angles will be enumerated twice, 
but the duplicates can be identified and removed during the sorting process. 

The usual sorting algorithms handle n elements in 0(nlog(?r)) operations, if we provide them a 
constant time comparison operator. The goal of this section is to explain how to compare two hinge 
angle via their generating triple coding. This comparison needs constant time. Using a quicksort on 
generated triples, the enumeration and the sorting of the hinge angles can therefore be computed 
in 0(m 3 log(m)). 

Lemma 4 (Integer-Based Comparison) Let a G Ah which is associated to the triple ( p , q, k ) 
and a' G Ah which is associated to the triple (p',q',k'), it is possible to decide, in constant time, 
whether a < ol using by knowing p, q, k,p ’, q’ , A/. 

Proof(sketch): With the notations previously introduced we have: cos(a) = (p(k + |) + 

q\)/(p 2 + q 2 ) and sin{a) = (pX - q(k + \))/(p 2 + q 2 ). 

The key-argument is that the numbers used in the expression of cos(a) and sin(a) are all integer, 
except A, but 4A is a square root of an integer, thus the other equations will simplify. 

Thus the signs of cos(a) and sin(a) can be easily computed, reducing to the problem to the 
sign of an integer inequality that involves only usual operations. Hence, the “quadrant” of a can 
be computed with integer arithmetic-only. 

If a and a' belongs to the same quadrant then we have to compare cos(a) and cos(a / ). This 
can also be reduced to determining the sign of an integer expression, that can be computed using 
only integer arithmetic. □ 


4.2 Moving from [r Q -] to [ r a +] 

Let a G An, we now know that exists z s (a ) G Z[i] and ^(^(a:)) + | G Z and also Q^z^a)) > 0. 
Due to the trigonometric orientation, and to the chosen discretization operator, lim [rv] ( z s ) = 

a'i— >a,a'<a 




[I’aK-^s) while lim [r a ](z s ) = ([?’ Q ](z s ) — l). However, we can also notice that lim [r Q /] (—z s ) 

Oi'\-^a.,a'>Oi cx'\—>ot,a'<a 

([ r a]( — z s) + 1) while lim [?v]( — z s ) = [r a ](—z s ). These results can be summarized in the fol- 

a'\—>a,a'>a 

lowing theorem: 

Proposition 4 MQ E {0,1,2,3}, lim [r Q+e ] ((z s i®) — [r a _ e ](z s i®)) = i® +2 

e\— >0,e>0 v ' 

Let if zo denote the map such that if z is in (^o^(2n+l)), with n, Q E N, then ip Z0 { z iP) = P+i® +1 
else z 0 (z,p ) = p. Thus, it can now be stated that for any z in Z[i],[r a +\{z) = ’fi Za { z i [ r a-](A)) 

4.3 Sketch of the Algorithm 

With the previous statements, we now have the necessary elements to understand the fundamental 
ideas of the algorithm. This version of the algorithm simply update a discretized rotation map and 
call a function that needs to be notified at each update. We will later see a more complex version 
which explains how to make rotate incrementally an image without storing any copy of the original. 
See Algorithm 1. 

Most of the subtilty of the algorithm actually proceeds from the mysteries of the sequence of 
hinge angles. (See Section 15.41) In real implementations, small technical details shall be solved: 
Naturally, some additional code has to be added to take in account the fact that most programs 
store images into rectangular buffers. Also, this version of algorithm actually jumps over the hinge 
angles, and it does not compute the p m (a) for the associated to the hinge angle it self. Of course, 
this is easily fixable: Apply ^ on the two first quarters, then call the procedure for notification; the 
algorithm has then to terminate the transformation associated to the hinge angle by computed the 

next two quarters, and to call once more the notification procedure. __ 

Algorithm 1: incremental_rotation_via_map_and_notifications(doextproc) 

1: An\m UstjandsortJiingejangles{rn) 

2: Vp E Z[i],R\p] <-p,(z 8 ,z d ) <- beg(A n \m(z s ,z d )) 

3: while true do 

4: (z s ,z d ) <- next_spc An \ m (z s , z d ) 

5: for Q <— 0 to 3 do 

6: k=0 

7: while |(2 k + l)z s | < m do 

8: op = R((2k + 1 )z s i®) 

9: np R((2k + l)z s i®) + 

10: R((2k + 1 )z s i®) <— np 

11: k i— k + 1 

12: end while 

13: end for 

14: doextproc(R, z s ) 

15: end while_ 

r denotes the radius of the image, op means old position, np stands for new position, the func¬ 
tion “doextproc” is a parameter function that is called each time the image has been set up in 
configuration that corresponds to the image of a discretized rotation. 

4.4 Application to Rotation of Images 

The previous version of the algorithm computes an image of the rotation map, however the dis¬ 
cretized rotation is not intrisically bijective on 7? and it is not suitable to compute incrementally 







rotations of an image without any copy of the original. 

It is well-known that a point of the discretized rotation has never more than two antecedents 
by discrete rotation. Thus, to create lossless discretized rotation, one natural idea is to store both 
antecedents when the function is not injective. This requires one additional layer: hence, an image 
in our algorithm shall be an element of Ql z W x I°> 1 }). 

Practically, the code is modified such that: the dataspace on which rotation are computed can 
support up to two “colors” for each position of Z[i]. We will use an arbitrary order (the lexical 
order or anything fast to compute) to decide which pixel will stand on the layer 1. 

The details that have been added allows to compute the 2-layers discretized rotation. See 
Algorithm 2. 

4.5 Analysis of the Complexity 

This last algorithm has really different complexity comparing to the usual algorithm for rotations. 
This algorithm is a bit slower than the usual one in the general case, but it is faster when incremental 
rotations are needed. Moreover this algorithm is slowed down by the fact for the moment we cannot 
assume that there exist a good structure to store the table of hinge angles. While, it is assumed 
that each processor uses a sine/cosine table which allows to compute sine and cosine very efficiently. 

Space Complexity : The second part of the algorithm uses the amount of memory required to 
store a 2 layer image, plus the rotation map. Thus, the algorithms uses about 3 K(m 2 ), (where K 
is the memory cost to store one color, or one vector (assuming that these two data-types can be 
stored with the same number of bytes K )). The traditional algorithm uses about 2 K(m 2 ) bytes 
of memory. Thus, for this part the space requirements are of the same order 0(m 2 ) and similar 
in terms of multiplicative factors. The first part of the algorithm requires 0(m 3 log(m)) bytes of 
memory to construct a b -tree of hinge angles in An\ m . This table can be computed once for all. 

Time complexity: The list-and-sort procedure for hinge angles requires 0(m 3 log(m)) operations. 
The time of one iteration of the loop in the second part is intrisically linear in m (the complexity 
of the user contributed function, doextproc can of course decrease these performances). Although, 
the main loop is called 0(m 3 ) times and contains another loop, the algorithm also requires only 
0(m 3 ) operations: we update only m 2 pixels and each pixel crosses at most 4m times the dual of 
the grid. To compute all p(a)\ m , for a £ [0, 2ir[ with traditional rotation algorithms would have 
needed 0(m 5 ) operations. 

Algorithm 2: incrementaljrotatiori-of-anjimage(imgi, doextproc) 


1 

An m list-andsortJhingejangles{m) 

15 

alone <— true; 

2 

Mp £ Z[i], R\p] <— p, L[p] <— 1 

16 

for d <— 0 to 3 do 

3 

(z s ,z d ) <- beg(An\m(z s , z d )) 

17 

if R\p + i d ] = op then 

4 

while true do 

18 

alone *— false; p' <— p + 

5 

(z s , z d ) «- nextjspc An \ m {z s , z d ) 

19 

end if 

6 

for Q <— 0 to 3 do 

20 

end for 

7 

k=0 

21 

if -i alone and L[j/] = 2 then 

8 

while |(2 k + l)z s | < m do 

22 

I- 'P'] I 

9 

p <— (2k + 1 )z s 

23 

suoap{img\ [op ], img 2 [op] ) 

10 

op <— R[p\ 

24 

end if 

11 

ol <— L[p\ 

25 

alone <— true; 

12 

np <— op + 

26 

for d <— 0 to 3 do 

13 

img 2 (np) «- img ol {op) 

27 

if R[p + i d ] = np then 

14 

R\p\ nv, L\p\ <- 2; 

28 

alone <— false; p' <— p + i 






29 

end if 


38 

else 

30 

end for 


39 

L[p\ «— 1 

31 

if -1 alone then 


40 

swap(img\ [np ], img -2 [np ]) 

32 

if p 1 < p then 


41 

end if 

33 

L\p'} v- 1 ;L[p] e 

- 2 

42 

end while 

34 

else 


43 

end for 

35 

L\p'} <- 2;L[p] - 

- 1 

44 

doextproc(R , z s ,img\) 

36 

swap(imgi [np], 

img 2 [np ]) 

45 

end while 

37 

end if 





The incremental discretized rotation algorithm. The notations are the similar to the one used 
for previous algorithm. Additionally, img\ and img 2 denotes the two layers of the picture. L is a 
part of the configuration that is used to memorize the destination layer of a pixel. 

4.6 Open-Source Implementation 

We have written an implementation of this algorithm. The C++ code can be downloaded from 
http: //perso. ens-lyon. fr/bertrand.nouvel/transitive-rotations/. This implementation 
relies on similar ideas but is actually slightly different: For historical reasons and other reasons 
that specific of the future of the code, it uses two types of configurations. The first kind is p a \m 
and is used to compute images by discrete rotation. While C a (see EEOE!) is used for some 
additional checking. 

5 Complement of specifications 

5.1 More on Angular Topology of Finite Balls in Z[i\ 

We denote by Iu\m the open intervals of angles delimited by two consecutive hinge angles of 
An\m- For any angle a G A and any m £ N either a G either there exists a unique 

I G In\m such that a € I. We define the function of Q m from A to A, such that for each hinge 
angle a G An\ m , Qmipt ) is the angle a itself; and for any angle a that belongs to an interval 
/ =]ao,ai[G In\m, we define Q m (a) as the Pythagorean angle which belongs to I which has the 
smallest radius. We define An\m = Qm(A). The set An\m has the following properties: Va G 
A, 3a' G A~}{ | m , such that p m (a') = p m (a) and Vai,a 2 G An\m, («i / a 2 ) =+ p m {of) / p m {pt 2 )- 
As p m can be seen as a bijection from An\m to p m (A), we define the function <f m that associates 
to each configuration p(a) the unique angle a! G An\m such that p(a) = p(a'). Note also that, if 
a, a' € I =]ao, aq[ and / n A r = 0 then p m (a) = p m (a) 

The hinge angles also have also more properties which are beyond the scope of this paper, and 
whose proof is mostly arithmetical: If a is a hinge angle then ka is a hinge angle for k 7 ^ 0 mod 3 
(by recurrence ka is a hinge angle then (k + 3)a is also a hinge angle) , if a is a non flat hinge 
angle then its sine and its cosine are rationally dependant, {Ze*“} forms span of 2 gcd(x,y) line of 
slopes | in (M/Z)[*]... 

Note: The set of hinge angle that has been obtained and studied here is dependant of the 
discretization that has been used (and of the real center of the rotation). 

5.2 Rotations of the Entire Grid Z 2 

If we consider infinite configurations, then Z[i] z M i s i n bijection with A. More precisely, for any 
a, a' € A, a = a', [r a \ = [r(J iff and a = ol . This is trivial since lim mrt0O (^m 0 Pm( a )) — a — i ► 0, . 









Moreover, there exists a convergent process that consists computing gradually each p m (a) with 
to incrementing through time from 1 to oo. This process is a well defined procedure to describe 
discrete rotations. 

5.3 e - Quasi- Transitivity 

The rotation algorithm specifies which map ^ should be applied to transform one rotation map 
[r a -](img) into another [r a +](img) (with respect to the transformations that have previously been 
applied). To create a transitive rotation, the cr-rotation-map of previously applied transformation 
needs to be embedded within the picture as a part of the data manipulated by the algorithm 1 . 

The incremental rotation has been implemented here as an endless process. The procedure 
rotates gradually a picture and regularly calls a function to notify the process that the image has 
a new layout. The algorithm can be easily modified to stop when it reaches a specified angle. 

From any rotation configuration C £ ((Z[?']) z [*H m ) it is possible to extract the associated slope 
to the configuration, by computing lim„ woo tan((j) m (p m (a))). From the embedded configuration 
C = p\m{oido) then we can restart the algorithm from a^ 0 = 4>m(C) and then we can rotate 
incrementally to a new angle a c ii • This process can be repeated without any incidence on the 
precision: the resulting map will stay a discretized rotation map. 

If we decide to denote by [r a ], the application of the previous algorithm from an angle ad 0 
(computed from the embedded configuration), to an angle a ^ = a.d 0 + a. Then, this operation 
would be e - quasi transitive, in the following meaning: |r a ] o |r^J = \r a+ p +e ^ a ^)\ . 

More general rules may be found later to apply rotations knowing the input but with broader 
constraints on the admissible type of permutations of Z [z]: We can imagine to allow rotations from 
various centers , translations... 

5.4 Open Question: Efficient Procedure for Computation of the Successor of 
an Hinge Angle 

One of the remaining problem of our algorithm is that the only way to compute efficiently the 
successor angle is now, up to our knowledge, to construct and use a table of the hinge angles. 

However, the question of the structure of hinge angles seems to be linked with some famous 
arithmetic and number theory results. Also, the hinge angles can be seen as the subset of the more 
general set Ms of angles such that a verifies: acos(a) + 6sin(a) = c, with a, b, c £ {—to, ..., to} 

For our algorithm, we are more interested by the structure induced on the generating triples 
( p , q , k ) than by the values of the angles themselves. Although the equation seems familiar due to 
its similarity to formulas for rotations, very few is known about its integer solutions. 

6 Conclusions and Perspectives 

This algorithm for discretized rotations has the numerous useful properties: Since it is an exact 
algorithm, that is valid on any size of datas, it returns the same result as the direct discretized 
rotation of for the associated angle. However, this algorithm proceeds incrementally, and it will be 
useful in all procedures that need to do some checking throughout the rotations. Additionally to 
be more exact, this algorithm is faster to compute than incremental rotations that would iterate 
the discretized rotation. More precisely it leads to 0(n 3 ) updates of pixel instead of 0(ji 5 ). One 

x Of course, it can also be imaginable to record the angle of which the image has been rotated before, but actually 
we think it is more applications will be interested by a configuration that describes the influence of previously applied 
transformations on the topology of the space 
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of its other property is to be very accurate in terms of spatial error since the error made on the 
position of an image is at most of (which is optimal). The fact that it does not use any sine 
nor any cosine strengthen this accuracy. 

This work seems to be extendable to the 3-shears rotations. The main advantage that motivates 
the use of similar methods to the 3-shears rotations is that these rotations are natively bijective. 
Therefore, the additional layer becomes obsolete; this would simplify the algorithm. Some refine¬ 
ments on its associated dynamical system associated with hinge angles are certainly required before 
to can initiate this study. 

Finally, although the hinge angles are sortable using only integer arithmetic, the procedure is 
slowed down by the use of large integer numbers. Thus, to speed up the algorithm, as suggested in 
Section IP1 a trend for ongoing research is to find an algorithm to return efficiently the successor 
of an hinge angle in An\m- 
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